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ABSTRACT 

A test program to evaluate a new transonic axial 
cortpressor stage was conducted. The stage was designed (by 
Nelson Sanger of NASA Lewis) relying heavily on CFD 
techniques while minimizing conventional enpirical design 
methods. The stage was installed in the NPS Transonic 
Coitpressor Test Rig and instrumented with fixed teitperature 
and pressure probes. A new PC-based data acquisition system 
was commissioned and programmed for stage performance 
measurements. These were obtained at 50, 60, 65, 70, and 80% 
of the design speed before failure of the spinner retaining 
bolt led to the loss of the stage. The flow through the 
rotor was analyzed and the rotor performance predicted using 
a 3-dimensional viscous code (RVC3D) . The predicted rotor 
performance agreed qualitatively and was nximerically 
consistent with the measured stage performance. 



V 



vi 






TABLE OF CONTENTS 






I . INTRODUCTION 

II. TRANSONIC COMPRESSOR RIG AND INSTRUMENTATION.. 



A. 

B. 



COMPRESSOR INSTALLATION..... 

I NSTRUMENTATI ON 

1. Measurement Devices.. 

2. Instrument Placement. 

TEST PROCEDURE 









III. DATA ACQUISITION AND REDUCTION 



A. 

B. 



C. 

D. 






OVERVIEW...... 

HARDWARE 

1. Scanivalve Control. 

2. Scanning Digital Voltmeter.. 

3. Counter/Totalizer., 

SOFTWARE 



DATA REDUCTION. 

IV. RESULTS OF THE TEST PROGRAM. 






A. INTRODUCTION. 

B. STAGE PERFORMANCE MAPS. 

C. PRESSURE PROFILES.. 

RESULTS OF CFD ANALYSIS. 






1 

7 

7 

8 

.8 

13 

14 






A. 

B. 

C. 

D. 



INTRODUCTION. 
GRID GENERATION. 
COMPUTATION. 
RESULTS.. 






VT. CONCLUSIONS AND RECOMMENDATIONS 
LIST OF REFERENCES 



.17 

.17 

.17 

.19 

22 

,.26 

..26 

.27 

,.31 

..31 

.31 

-.33 



— 47 

.47 
.47 
.49 
50 

55 

57 



vii 



APPENDIX A 



ANALYSIS OF THE SPINNER BOLT FAILURE. 



59 



A. 

B. 

C. 



OVERVIEW. 

CONSTRUCTION. 

ANALYSIS 






59 

59 

60 



APPENDIX B 



RESULTS.. 









..67 



APPENDIX C SOFTWARE USER'S GUIDE 



.71 



A. 

B. 



HPVEE., 

RVC3D., 









71 

98 



INITIAL DISTRIBUTION 



LIST. 



.105 



viii 



LIST OF FIGURES 

1. Schematic of the Transonic Compressor Test Rig 3 

2 . View of the Assembled Stage 4 

3 . View of the Instrumented Assembly 

With Lucite Case Wall 5 

4 . Schematic of the Flow Path 

Showing Measurement Stations 10 

5 . Torque Calibration Curve 12 

6. Schematic of the Data Acquisition System. 18 

7. Schematic of Switchbox Connections to HG78 20 

8. Schematic of Differential Tenperature Measurement 25 

9. Pressure Ratio (T-T) vs Referred Mass Flow Rate 35 

10. Adiabatic Efficiency (Tjj) vs Referred Mass Flow Rate 36 

11. Adiabatic Efficiency (TiJvs Referred Mass Flow Rate 37 

12. Adiabatic Efficiency (T|j) vs Referred Mass Flow Rate 38 

13. Referred Torque vs Referred Mass Flow Rate 39 

14. Stage Temperature Rise vs Referred Mass Flow Rate 40 

15. Ref erred Horsepower vs Referred Mass Flow Rate 41 

16. Exit Pressure Distribution at 70% Speed 42 

17. Exit Pressure Distribution at 80% Speed- 43 

18. Axial Static Pressure Distribution at 70% Speed 44 

19. Axial Static Pressure Profile at 80% Speed. 45 

20 . Computational Grid for RVC3D Calculations. 48 

21. Rotor Total Pressure Ratio vs Corrected (Referred) Flow 

Rate Computed Using RVC3D.... 52 

22 . Rotor Adiabatic Efficiency vs Corrected 

(Referred) Flow Rate Computed Using RVC3D 53 



IX 



1 



X 



LIST OF TABLES 



1. Design Parameters for the Sanger Stage 2 

2 . Location of Probes at the Stage Exit 13 

3. Scanivalve Port Assignments 21 

4 . Scanning Multiplexer Channel Assignments 23 



XI 



ACKNOWLEDGEMENT 



The author would like to acknowledge the financial 
support of NAWCAD, Propulsion and Power Engineering, Code 
4. 4. 7.1, for the experimental study reported here. The work 
was performed under Contract N004219WR01221 . 



I . INTRODUCTION 



A new transonic axial compressor stage designed 
specifically for the Naval Postgraduate School (NPS) 
Transonic Coitpressor Rig (TCR) was conpleted at NASA Lewis 
in 1994, by Nelson L. Sanger [Ref. 1]. A schematic of the 
TCR is shown in Fig. 1. The new stage design was developed 
to serve as a vehicle for instructional purposes, and to 
provide a test bed for research. With the emergence of 
Conputational Fluid Dynamics (CFD) as a reliable tool, the 
design of the stage developed into, as Sanger described it, 
"...the ultimate CFD validation experiment" . The design 
parameters produced from this process are presented in Table 
1 . 

The manufactured stage, after assembly in the TCR, is 
shown in Fig. 2. The stage was installed in the same manner 
as previously operated stages and as documented in Ref. 2. 
However, a Lucite case wall was installed, as seen in Fig. 
3, and a new data acquisition system was employed. The data 
acquisition system used a Personal Conputer (PC) riinning 
Hewlett-Packard HPVEE software to control various 
instruments in a HP75000 VXI-bus mainframe. The stage was 
tested and measurements taken at' 50, 60, 65, 70 and 80 
percent of design speed. Testing was vinavoidably halted 
during the test at 80% speed due to the failure of the 
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spinner retaining bolt, which resulted in the loss of the 
stage . 

The present report docxaments the first test program of 
the Sanger stage. The report contains three main sections. 
The first section details the rig and the newly developed 
data acquisition system. The second section focuses on the 
results . The final section gives the results of a CFD 
analysis of the rotor using the "Rotor Viscous Code 3 
Dimensional" (RVC3D) code developed by R. Chima [Ref. 3]. 



Parameter 


Final 

Selection 


Rotor Pressure Ratio 


1.61 


Stage Pressure Ratio 


1.56 


Tip Speed 


1300 ft/sec 


Design Mass Flow 


17.09 Ibm/ sec 


Specific Mass Flow 


35 Ibm/ sec -ft^ 


Specific Head Rise 


0.246 


Tip Inlet Relative Mach 
Number 


1.28 


Hub/ Tip Radius Ratio 


0.51 


Rotor Inlet Ramp Angle 


28.2 degrees 


Power Required 


457 HP 



Table 1. Design Parameters for the Sanger Stage. 
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Fig. 1 Schematic of the Transonic Coirpressor Test Rig. 
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Fig. 2 View of the Assembled Stage. 
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Fig. 3 View of the Instr\amented Assembly With Lucite Case 

Wall . 
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« 
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II. TRANSONIC COMPRESSOR RIG AND INSTRUMENTATION 



A. COMPRESSOR INSTALLATION 

The general layout of the TCR is shovm in Fig. 1. The 



manner in 


which 


the 


stage was 


mounted 


on 


the rig was 


unchanged 


from 


the 


previously 


installed 


Vavra stage 


docxamented 


in Ref - 


. 2. 


The flow entered the 


rotor axially. A 



honeycomb section ensured that the flow left the conpressor 
with zero circumferential velocity. In this manner the 
torque supplied to the rotor was equal to the torque 
experienced by the blading and honeycomb of the stator 
section. The stator section was mounted on ball races and 
was free to rotate against a flexure instrumented with 
strain gauges. The strain gauge reading was then a measure 
of the torque supplied to the rotor. 

The tip clearance between the rotor and the steel 
casing that had been used previously was a concern. The 
decision was made to use a similar clear plastic (Lucite) 
casing xantil the steel casing could be modified to have an 
abradable material where the rotor might rub. The Lucite 
casing proved to function extremely well. The static tip gap 
clearance was 0.010-0.012 inches which was larger than the 
0.004 inches the design called for; however, this ensured 
that no rxib would occur during the first test of this unique 
stage. Additionally, the transparent case provided an 
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opportunity to see transonic turbomachinery in operation. 
One interesting observation was that excess oil from the 
bearing lubrication system indicated a flow pattern and was 
seen to collect between the casing and the rotor tip. 

B. INSTRUMENTATION 

1 . Measurement Devices 

Fig. 4 shows the location of measurement stations 
relative to the stage. The instrumentation for stage 
performance measurements included Kiel pressure probes, 
Kiel /thermocouple combination probes, static pressure taps, 
torque flexure, flow nozzle, and a magnetic pick up. Total 
pressure was measured using United Sensor Kiel Probes. Two 
types were used, a 1/16 inch "Miniature Head" Kiel Probe 
(United Sensor KAA-8, Ref. 4) and 1/8 inch "Standard Head" 
combination Kiel/ thermocouple Probe (United Sensor KT-8-J- 
12-C, Ref. 4) . The 1/16 inch mini probe had a yaw range of 
52° and pitch range of 40°. Range was defined by the 
manufacturer as the point at which error was equal to 1% of 
dynamic pressure. The 1/8 inch standard probe had a yaw 
range of 48° and a pitch range 45°. The static pressure ports 
were 1/32 inch taps located in the wall along the length of 
the casing, and in the hub of the stator at stations 2 and 
3. 
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Thermocouples were J type (iron-cons tantan) in a probe 
combined with a Kiel type impact probe. The thermocouple had 
an exposed junction and bleed slot. 
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Fig. 4 Schematic of the Flow Path Showing Measurement 

Stations . 
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The torque acting on the stator and exit-honeycomb 
section (equal to the rotor torque) was measured by the 
deflection of a cantilevered flexure. The cantilever arm was 
instriimented with four strain gauges providing a near linear 
response to bending. The arrangement and calibration method 
of the torque measurement system is extensively documented 
in Ref. 2. Calibration of the torque measuring system was 
necessary each time the rig was rebuilt. The calibration 
curve obtained in the present study and resulting equation 
used in the data reduction program are shown in Fig. 5. 

Rotational speed was measured using a magnetic flux 
cutter and a disk with six holes mounted on the shaft of the 
drive turbine (Fig. 1) . This device produced six pulses per 
revolution making the frequency (in Hz) from the pickup 
multiplied by 10 equal to the revolutions per minute (RPM) of 
the stage. 

Mass flow rate was measured using a flow nozzle located 
upstream of the stage as shown in Fig. 1. Mass flow rate was 
calculated using the pressure difference of the total 
pressure (p^) and static pressure (p^) at the flow nozzle and 
the total temperature of the flow. The differential pressure 
between p^ and p^ was displayed on a water manometer in the 
control room to allow precise control of mass flow rate 
settings during the throttling process. 
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Fig. 5 Torque Calibration Curve. 
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2. Instrument Placement 

Two combination probes were placed upstream of the 
rotor, each two inches from the wall of the casing to 
measure inlet total pressure and temperature. At station 3, 
17 Kiel probes and three Kiel/ thermocouple combination 
probes were placed as shown in Table 2 . 



Probe 

# 


Probe 

Type 


Radial Distance 

(y) 


Angle 




18 


Kiel 


0.7 


15 


15 


Kiel 


0.4 


30 


12 


Kiel 


0.1 


45 


41 


Static Port 


n/a 


60 


11 


Kiel 


0.05 


75 


29 


Kiel 


1.8 


90 


/ 


Positive \ 

Angle \ 


26 


Kiel 


1.5 


105 


23 


Combo 


1.2 


120 


20 


Kiel 


0.9 


135 


17 


Kiel 


0.6 


150 






14 


Kiel 


0.3 


165 


+y inches 


42 


Static Port 


n/a 


180 




plug 


n/a 


n/a 


195 




28 


Kiel 


1.7 


210 


25 


Kiel 


1.4 


225 






22 


Combo 


1.1 


240 




19 


Kiel 


0.8 


255 


View LOOKing uownsueam 


16 


Kiel 


0.5 


270 


13 


Kiel 


0.2 


285 


43 


Static Port 


n/a 


300 


30 


Kiel 


1.85 


315 


27 


Kiel 


1.6 


330 


24 


Kiel 


1.3 


345 


21 


Combo 


1 


360 



Table 2. Location of Probes at the Stage Exit. 



13 



The pressures acquired from the impact probes at 
station 3 (probe numbers 11-30 in Table 2) were mass 
averaged to obtain the total pressure at the exit of the 
stage. 

Static pressure ports were placed at stations Inf, 1, 2 
and 3 along the Lucite casing (Fig. 4). Additionally, static 
taps were drilled in the stator to measure hub static 
pressure at station 2 . The exit spacer between the stator 
and honeycomb also had four static taps to measure exit 
static pressure on the hub at station 3 . 

C . TEST PROCEDURE 

The test compressor was driven hy two opposed-rotor air 
turbine stages (on a common shaft) supplied from the 
laboratory cortpressed-air power supply. The laboratory 
conpressed-air power supply was a 12-stage Allis-Chalmers 
axial compressor supplying up to 11 pounds per second at up 
to two atmospheres gauge pressure. The test corrpressor was 
connected by a sp lined shaft to the rotor of the turbine 
drive lanit. The speed of the drive unit was controlled using 
an electrically actuated butterfly valve. Air was drawn into 
the test compressor through a rotating plate throttle 
followed by a settling chamber containing perforated plates 
(Fig. 1) . A flow nozzle was positioned upstream of the 
rotor. The calibration of the flow nozzle was documented in 
Ref. 5. 
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Preparation for a run consisted of bringing the Allis- 
Chalmers cortpressor up to stable speed, verifying that 
compressed air was available to the balance piston (which 
maintained low axial loading on the compressor bearings) , 
and setting the lubrication system. The bearing liibrication 
system was an oil-mist system. The operating temperature of 
the bearings was found to be directly related to the supply 
pressure of the air driving the oil mist. Increasing the 
pressure of the oil mist supply air had the observed effect 
of lowering the bearing operating temperature. Adjusting the 
oil droplet rate down to 2-3 drops-per-minute was essential. 
According to the bearing manufacturer, wetting the bearing 
well before a run was of extreme inportance. Allowing the 
oil mist to operate for fifteen minutes before operating was 
the procedure adopted. 

During a conpressor test, operating data were monitored 
by the operator at the control console. Records were kept of 
oil drop rate and pressure to the bearing lubrication 
system, of the tenperatures of the turbine and conpressor 
bearings, and of vibration levels from four accelerometers 
on the turbine drive unit. The axial force was constantly 
monitored and the balance-piston air adjusted to keep the 
axial force constant at a level approximately equal to the 
force produced at 5000 RPM. These data were mainly of 
interest for the safety and operation of the unit and will 
not be reported. 
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A typical test would be conducted at a constant speed. 
The run started with the compressor throttle wide open. 
After the speed had settled to the desired level, a 
measurement recording cycle was started. Then the throttle 
would be closed fractionally until the flow was reduced to 
the desired amount. Then the drive turbine speed was 
adjusted until it was steady again and at the desired level. 
This procedure was repeated until stall occurred. After the 
stall point was reached (up to 70% speed only) the process 
was progressively reversed to obtain intermediate settings 
until the throttle was fully open. 
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III. DATA ACQUISITION AND REDUCTION 



A. OVERVIEW 

In order to standardize equipment for data acquisition 
in the three buildings at the NPS Turbopropulsion Laboratory 
the decision was made to use the HP75000 Series B VXI-Bus 
Mainframe controlled by HPVEE software running on a personal 
computer (PC) . This allowed the Hewlett-Packard (HP) 
instmments with the HP-IB interface that had been used 
previously to remain viable while incorporating a new state- 
of-the-art data acquisition system. In the present work the 
mainframe with an internal digital voltmeter (DVM) , two 
scanning multiplexers, a switchbox multiplexer, a counter, 
and a Quad 8-bit Digital I/O Module were used exclusively. 
The system, shown in Fig. 6, proved to be flexible, easy to 
interrogate when connecting new sensors, and provided real 
time feedback when acquiring data. 

B. HARDWARE 

The use of the HP75000 Mainframe allowed a variety of 
instruments to be grouped together and addressed directly. 
Communication between the PC (with a HP 82341C Controller 
Card installed) and the mainframe was via a HP-IB (IEEE-488) 
interface cable. The functions performed with the 
instruments in the mainframe were Scani valve control, 
measurement of Scanivalve transducer output, thermocouple 
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measurement , 
counting . 



strain gauge 



measurement 



and 



frequency 




Fig. 6 Schematic of the Data Acquisition System. 
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1 . 



Scani valve Control 



The Scanivalve and Scanivalve control (HG-78) system 
that the present system was designed to replace, are 
described in Ref. 6. Scanivalve control involved stepping 
and homing the 48-port pneumatic valve and ensuring that the 
correct port was selected and measured. One HP E1345A 16- 
channel multiplexer was configured as a switchbox as shown 
in Fig. 7. This configuration allowed the control of four 
different Scanivalves with one switchbox. The measurements 
required for the test program required the control of only 
one Scanivalve at position 1. To reset the Scanivalve, 
channel 12 was closed, sending the signal to the HG-78 to 
reset the Scanivalve connected to position 1. To step to the 
next port, channel 8 was closed, sending the signal to 
advance one port. The position of the Scanivalve was 
obtained from the HG-78, which contained LED displays of the 
port setting, manual step and home controls (in parallel), 
and signal conditioning for the Scanivalve transducers. From 
the HG-78, eight TTL lines (from each Scanivalve) were 
connected to one of four addressable sections of the HP 
E1330B Quad 8-bit Digital I/O card. The port which the 
Scanivalve was currently set at was then read (using 
negative logic) and converted to a decimal nixmber using 
Binary Convert Decimal (BCD) logic. The assignment of the 48 
Scanivalve ports is shown in Table 3 . 
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voltage sense 

TREE 

VOLTAGE SENSE 
DIRECT 



6 PIN CANNON 




Fig. 7 Schematic of Switchbox Connections to HG-78. 
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Port# 


Scanivalve 1 


1 


Tare 


2 


Scale 


3 


Flow Nozzle P6 


4 


Flow Nozzle Ps 


5 


Inlet Ptinfinily 1 Q1 3 O'clock 


6 


Inlet Ptinfinity 2 Q1 9 O'clock 


7 


Pnozzle static 


8 




9 




10 




11 


Kiel Probe 


12 


Kiel Probe 


13 


Kiel Probe 


14 


Kiel Probe 


15 


Kiel Probe 


16 


Kiel Probe 


17 


Kiel Probe 


18 


Kiel Probe 


19 


Kiel Probe 


20 


Kiel Probe 


21 


Kiel/TC Combination 


22 


Kiel/TC Combination 


23 


Kiel/TC Combination 


24 


Kiel Probe 


25 


Kiel Probe 


26 


Kiel Probe 


27 


Kiel Probe 


28 


Kiel Probe 


29 


Kiel Probe 


30 


Kiel Probe 


31 


Hub Static Pressures P2 1 


32 


Hub Static Pressures P2 2 


33 


Hub Static Pressures P2 3 


34 


Hub Static- Pressures P2 4 


35 


Hub Static Pressures P3 1 


36 


Hub Static Pressures P3 2 


37 


Hub Static Pressures P3 3 


38 


Hub Static Pressures P3 4 


39 


Shroud Static Pressures PI 3 O'clock 


40 


Shroud Static Pressures PI 9 O'clock 


41 


Shroud Static Pressures P3 3 O'clock 


42 


Shroud Static Pressures P3 6 O'clock 


43 


Shroud Static Pressures P3 9 O'clock 


44 


Shroud Static Pressures Old style 1 P3 


45 


Shroud Static Pressures Old style 2 P2 


46 


Shroud Static Pressures Old style 3 PI 


47 


Shroud Static Pressures Old style 4 PI 


48 


Shroud Static Pressures Old style 5 Pi 




Time References are from turbine out 



Table 3 . Scan! valve Port Assignments 
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The logic sequence for acquiring data from the 48 ports 
of the Scani valve was as follows : 

1. Reset to port 1 

2 . Ensure the Scanivalve was actually at port 1 

3 . Take the transducer reading using the DVM 

4 . Step to the next port 

5 . Ensure the Scanivalve was actually at the next port 

6. Take the transducer reading using the DVM 

7 . Go to step 4 until all 48 ports were read 

2. Scanning Digital Voltmeter 

Connecting a 16-channel multiplexer to the DVM using 
the 'Analog Bus' allowed the multiplexer and DVM to operate 
together as a Scanning Multimeter. ['Scanning Multimeter' 
implies that one or more multiplexers were used with the 
multimeter as part of the same instrioment using the same HP- 
IB secondary address] . The power of this configuration was 
that the multimeter was extended from having one set of 
terminals to having sixteen for each connected multiplexer. 
Two types of scanning multiplexer were eitployed. The first 
type was a thermocouple relay multiplexer module (HP 1347A) . 
A thermistor was contained inside the terminal module, and 
its resistance used to determine tenperature with which to 
compensate for temperature measurements made with the 
thermocouples . The second type was a standard relay 
multiplexer module (HP 1345A) , which allowed the measurement 
of sixteen different channels using a single multiplexer/DVM 
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connection. The channel assignment for the two multiplexers 
is given in Table 4. 



Multiplexer 

Channel 


Thermocouples /Mux 100 


DVM Scanner /Mux 200 


00 




Scanivalve 1 
Transducer 


01 


T^ Station Inf @ 3 
o ' clock 


Scanivalve 2 
Transducer 


02 


T^ Station Inf @ 9 
o'clock 


Scanivalve 3 
Transducer 


03 


T^ Station 3 Probe 21 


Scanivalve 4 
Transducer 


04 


T^ Station 3 Probe 22 


Torque-Flexure Strain 
Gauge Array 


05 


T. Station 3 Probe 23 




06 






07 






08 






09 






10 


At^ - Station Inf @ 3 
o'clock and Probe 21 




11 


AT^ - Station Inf @ 3 
o'clock and Probe 23 




12 






13 






14 






15 







Table 4. Scanning Multiplexer Channel Assignments. 

Readings from the four controlled Scanivalve 
transducers and the torque- flexure strain gauges were read 
directly and controlled by the DVM driver software. 

To measure temperature directly, the thermocouple wires 
were connected directly to the thermocouple terminal module 
and the voltage measured and converted to degrees 
Fahrenheit. To measure the instantaneous differential 
teirperature between stations 3 and Inf, a different process 
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was required. The thermocouples were connected as shown in 
Fig. 8. 

The logic process was as follows ; 

1. Measure the temperature at station Inf (using the 
thermocouple multiplexer) . 

2 . Convert the tenperature measured in step 1 to volts 
using the equation relating temperature to voltage for a J 
type thermocouple . 

3 . Measured the voltage between the two thermocouples 
(using the scanning multiplexer) . 

4 . Added the measured voltage of step 3 to the voltage 
obtained in step 2 . 

5. Converted the voltage obtained in step 2 to a 
tenperature using the equation relating voltage to 
temperature for a J type thermocouple. 
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TC Multiplexer TC Multiplexer 




IRON CONSTANTAN 
THERMOCOUPLES 



Fig. 8 Schematic of Differential Tertperature Measurement. 
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3. 



Coimter/Totalizer 



The Counter/Totalizer was used to measure the frequency 
produced by the magnetic speed pickup. The pickup produced 
six pulses per revolution. Therefore, the frequency 
measured, in Hz, multiplied by 10 produced the RPM of the 
rotor. The signal produced was extremely noisy near the zero 
volt level. To reject this noise, the level to trigger the 
counter was set to +1 volt. The rotational velocity of the 
rotor was verified using a calibrated strobe light. 

C. SOFTWARE 

HPVEE was the software used to control the 
instrumentation. Software drivers were available in HPVEE 
for all instruments inside the mainframe. The 'panel' driver 
for the specific instrument served two pu3rposes. First, it 
allowed the definition of a measurement state, that is the 
settings of all instrument functions. When the panel driver 
operated, the corresponding physical instrument was 
automatically programmed to match the settings defined in 
the panel driver. Secondly, the panel driver acted as a 
virtual control panel for interactively controlling the 
instrument . 

In addition to controlling the instruments, HPVEE also 
allowed immediate reduction of the measured values to 
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engineering units. This allowed immediate feedback of data 
and made troubleshooting simpler. 



D. DATA REDUCTION 



The data were reduced using HPVEE. In the description 
which follows, the final form of the data reduction program 
is described. The results presented were calculated with the 
following methods. 

The mass flow rate (in Ibm/sec) through the flow 
nozzle, was given by 



I P, A 

m = 33l35C-J—^—^- 
V 459.7 + 7 ,,^ 



1-0.021926- 



( 1 ) 



where A^=(p—P^) was in inches of water and the P. was in 

inches of mercury. The constant C represented the discharge 
coefficient and was calibrated to be 1.03. 

The first measure of the power absorbed by the 
compressor was obtained from the measurements of torque, M, 
and speed, N, as 

//Pci =1.5865 X 10'^ -M -A( (2). 

The second measure of power was obtained from the flow 
rate and teirperature rise across the stage, namely 

//Pc2 = 1.4137 3 -rj (3). 

The performance of the compressor was described in 
terms of referred quantities that retain their original 
units : 
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. # 






5 





( 4 ) 




where 6 = - 7 ^ 

Tarr 



and 5 = — 



and where T^^£=518.7°R and 



REF 



P^^j=29.92 inches of mercury. 

The performance to be evaluated was of the rotor and 
stator combination. The honeycomb flow straightener in the 
outlet duct was necessary in order for the torque to be 
measured, but the losses occurring in this section are not 
included in appraising the compressor performance. The 
total-to-total (T-T) pressure ratio and efficiency from 
station 1 to station 3 were the parameters which gave the 
performance of the stage . 

The total-to-total efficiency was defined as the ratio 
of the power required in an isentropic compression to the 
same outlet pressure, to the actual power required. The 
actual power was given by equation ( 2 ), using measurements 
of torque and speed, or by equation (3), using measurements 
of mass flow rate and temperature rise. The ideal power was 
given by 
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////>, =1.4137mC,((7;j)_-r,„,] (5) 

where is related isentropically to the pressures by 

T P 

^ \ /inf 

where y is the ratio of specific heats. Substitution into 
equation ( 5 ) , gives 




IHP,=U\31mC^T,^, 




(7) . 



Three different methods were used to calculate the 
total-to-total efficiency. Method 1 simply used the 
temperature rise and total-to-total pressure ratio across 
the stage 




( 8 ) 



The second method used inlet static pressure 
measurements and mass flow rate to calculate an inlet 
stagnation pressure. Defining 






P:. 



r-1. 



inf 






m 



Area 



(9) 



where the area corresponded to where the measurement was 
taken, it follows that 
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( 10 ) 




and Kinf ^p '^tinf • ^inf was given by 



A:„ = Viy^-s, 



( 11 ) . 



Thus 




( 12 ) 



was an efficiency which involved the mass flow rate 
measurement, whereas r|i did not. 

The final method for calculating efficiency depended on 
the measurement of torque and mass flow rate but did not 
depend on the measurement of the tenperature rise across the 
stage , namely 




( 13 ) . 
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IV. RESULTS OF THE TEST PROGRAM 



A. INTRODUCTION 

The performance of the stage was measured at 50%, 60%, 
65%, 70%, and 80% of the design speed of 27,085 rpm. The 
instrumentation remained constant during all runs and no 
changes to the data acquisition or reduction system were 
required. Some changes were made in how the reduced data 
were saved making plots of torque, horsepower and stage 
ten^erature rise versus mass flow rate only available for 
70% speed and 80% speed. Each set of plotted experimental 
performance data was fitted with a second order least- 
squares curve. The exit impact pressure profiles and axial 
static pressure distributions were also plotted. 

B. STAGE PERFORMANCE MAPS 

The total- to-total pressure ratio, efficiency 
calculated by equations (8), (12) and (13), referred torque, 
referred horsepower and stage temperature rise were plotted 
versus the referred mass flow rate, for constant speed 
conditions expressed as a percentage of design speed. 

Figure 9 shows the total- to- total pressure ratio vs 
referred mass flow rate. The plot shows the characteristic 
increase of pressure ratio and decrease of mass flow rate as 
the flow was throttled. At all speeds except for 80% speed 
the stage was throttled to stall. The loss of the stage 
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occurred before stall was reached at 80% speed. The open- 
throttle 'choke' condition for the stage could not be 
reached due to the pressure drop across the honeycomb flow 
straightener . 

Figures 10, 11, and 12 are efficiency plots. The 
ef f iciency T|, was expected to be the most accurate as it used 
torque and speed measurements which provided an integral 
measurement of the work to the compressor. The other two 
efficiency measurements used the temperature rise across the 
stage. Since the total temperature was not expected to be 
constant radially at the exit and the measurement was taken 
at only three points in the center of the passage, the r|, 

and Tlj measurements of efficiency were likely to be 
optimistic. The results show that rij varied in the same way 
as the other two efficiency measurements, but was 
consistently lower in magnitude. The data for rjj at 80% 

speed measured 2.5% lower than the other constant speed 
lines. This was possibly the result of drift in the strain- 
gauge bridge output (caused by increased temperature) which 
was not balanced to zero before the final run was started. 
The drift was observed after the loss of the stage. The 
results for and support this contention since no 
appreciable drop was observed in the 80% speed results. 
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The CFD analysis results shown on the performance maps 
were obtained for the rotor geometry, and did not include 
the stator. These results are discussed in the following 
section. 

Figure 13 shows the increased torque required as the 
compressor was throttled. Figure 14 shows the increase in 
the stage temperature rise as the compressor was throttled. 
The stage temperature rise was calculated as the average of 
the differential temperature of one probe at station Inf 
with two probes at station 3 . Figure 15 shows a comparison 
of the horsepower required to drive the stage calculated 
using the mass flow rate and stage temperature rise, (Eq. 
(3)), with that calculated using the torque and speed 
measurements (Eq. (2)). The horsepower measured by torque 
and speed was consistently greater than horsepower deduced 
from mass flow rate and temperature rise. 

C. PRESSURE PROFILES 

Figures 16 and 17 show the exit total pressure profiles 
at 70% speed and 80% speed respectively. (Note that the 
probes were at different peripheral settings in the radial 
direction shown in Table 2). The ratio of each measured 
total pressure at station 3 to the total pressure at station 
1 was calculated. Three different throttle settings are 
shown at each speed. At design conditions, the stage was 
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designed to have nearly uniform total pressure across the 
exit [Ref. 1] . The plots show that as the stage was 
throttled, at an intermediate throttle setting the total 
pressure was nearly uniform across the annulus, except near 
the walls. The plots also show that as stall was approached, 
which corresponded to the lowest mass flow rate curve, a 
separation had progressed radially outward from the hub. 

Figures 18 and 19 are plots of the static pressure 
measured at different axial stations at the same rotor 
speeds and throttle settings shown in the data in Figs . 16 
and 17 respectively. All pressures are scaled by the 
pressure at station infinity. 

At station 2 a strong static pressure gradient existed 
from hub-to-tip at all measured states of throttling and 
wheel speed. This gradient is consistent with the high 
tangential velocities expected at the rotor exit. The 
pressure rise through the stator behaved in a consistent way 
with throttling at each of the speeds measured. First, the 
pressure rise at the hub was large at all throttle settings. 
At the tip, the rise improved with closing throttle, until 
the static pressure difference between hub and tip at 
station 3 disappeared as the stage approached stall. 
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Fig. 9 Pressure Ratio (T-T) vs Referred Mass Flow Rate. 
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Fig. 10 Adiabatic Efficiency (Hj) vs Referred Mass Flow 

Rate . 
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Fig. 11 Adiabatic Efficiency {r\^) vs Referred Mass Flow 

Rate . 
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Fig. 12 Adiabatic Efficiency vs Referred Mass Flow 

Rate . 
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Fig. 13 Referred Torque vs Referred Mass Flow Rate. 
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Fig. 
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14 Stage Teit^jerature Rise vs Referred Mass Flow Rate. 
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Fig. 15 Referred Horsepower vs Referred Mass Flow Rate. 
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Fig. 16 Exit Pressure Distribution at 70% Speed. 
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Fig. 17 Exit Pressure Distribution at 80% Speed. 
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Fig. 18 Axial Static Pressure Distribution at 70% Speed. 
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Fig. 19 Axial Static Pressure Distribution at 80% Speed. 
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V. RESULTS OF CFD ANALYSIS 



A. INTRODUCTION 

Before starting the experimental work, the flow through 
the rotor was analyzed and the rotor performance was 
predicted using CFD. Specifically, Rotor Viscous Code 3-D 
(RVC3D) , a computer code for analysis of three-dimensional 
viscous flows in turbomachinery [Ref. 3], was used to 
compute performance maps of the rotor from 50% to 100% of 
the design speed. Additionally, using the graphics package 
FAST from NASA Ames, a number of parameters were calculated 
and graphically displayed for 100% speed. The predicted 
results for the rotor were compared to the experimentally- 
obtained maps for the stage and they were observed to be 
numerically consistent with each other. 

B. GRID GENERATION 

The geometry data for the hub, tip and blade for the 
stage were provided electronically by NASA Lewis. The data 
were in MERIDL (z, r, theta, dtheta) format. 

The geometry file was entered into Turbomachinery C 
GRID (TCGRID) which generated a three-dimensional C type 
grid conpatible with the Rotor Viscous Code 3-D (RVC3D) 
[Ref. 7]. The grid which was 185X40X49 (in i,j, and k) is 
shown in Fig. 20. The output grid was stored in plot3d 
multiple grid format. A FORTRAN program was used to rotate 
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the grid one grid width to form the multi-block grid that 
was an input to FAST. 




Fig. 20 Computational Grid for RVC3D Calculations 



48 




C . COMPUTATION 

The code solved the thin-layer Navier-Stokes equations 
with an explicit four stage Runge-Kutta scheme. The number 
of stages was optional; however, four was recommended by the 
author. A spatially-varying time step (At) was used to 
accelerate the convergence of the code. At was updated 
every fifty iterations. 

Since the flow was supersonic, a second-dif f erence 
artificial viscosity term was used. The value was set to 1.0 
in order to produce a well defined shock with little 
smearing or waviness. No first-order artificial viscosity 
term was required. 

RVC3D used an implicit residual smoothing routine which 
allows the multistage scheme to run at a high, unstable 
Courant number, but stability was maintained by smoothing 
the residual using an implicit filter. The initial Courant 
number chosen was 5.0. 

The code ran conpletely non-dimens ionally, with 
variables referenced to speed of sound, reference 
temperature and reference pressure. The input parameters 
that controlled the flow were the normalized blade row 
rotational speed per unit length and the ratio of the hub 
exit static pressure to the reference total pressure. 
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D. RESULTS 

Version 208 of the code was rian on three different 
machines, a Silicon Graphics Indigo2 workstation, the NPS 
Cray, and NASA Ames Cray. On the NASA Ames Cray, each 
iteration took 2.3 seconds of CPU time and required 9Mw of 
storage. The code was conpiled using FORTRAN 90 and 02 
optimization that gave good scalar and vector optimization 
and performed autotasking. Six processors were normally 
used. In order for the grid, solution, and restart files to 
be compatible with each machine, IEEE format was selected 
for the data input and output . 

To compute the operating lines at a constant rotational 
speed, from the choke region to the stall region, the ratio 
of the hub exit static pressure to the reference total 
pressure was increased in steps. For each p^^, at which 
a solution converged, the mass flow rate, the mass-averaged 
inlet and exit total pressures, and the adiabatic efficiency 
were recorded from the output file for plotting. When a 
converged solution could not be obtained {near stall) , the 
rotational speed was changed and a new range of values of 
Prat selected and entered. The operating speeds computed 
in this manner were 70%, 80%, 90% and 100% of design speed. 
Additionally, one point was computed at 50% of design speed. 
The results are shown plotted in Figs. 21 and 22. The 
results were observed to be extremely close to the results 
Sanger obtained using an Euler solver with bo\andary layer 
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modeling (Denton's Code [Ref. 1]). At design mass flow rate 
of 17.1 Ibm/sec Sanger obtained a total pressure ratio of 
1.61, and using RVC3D a total pressure ratio of 1.601 was 
found. 

The predicted results for the rotor at 70% and 80% 
speed are shown plotted with the experimental data for the 
stage in Figs 9 and 10 . A direct nvtmerical comparison is not 
meaningful since the predicted results were only for the 
rotor. However, the computed results are qualitatively 
consistent with what was measured for the stage behavior. 
The pressure ratio computed for the rotor was consistently 
higher than was measured for the stage, and the range of 
mass flow rate predicted at each rotor speed was consistent 
with the stage measurements. At 80% speed, solutions were 
obtained closer to choke (more open throttle) than could be 
achieved in the test program for the stage. 
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Fig. 21 Rotor Total Pressure Ratio vs Corrected (Referred) 
Mass Flow Rate Computed Using RVC3D 
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Fig. 22 Rotor Adiabatic Efficiency vs Corrected (Referred) 
Mass Flow Rate Computed Using RVC3D 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A test program was conducted to 80% of design speed to 
evaluate the performance of a new transonic axial coirpressor 
stage. Tests to 100% speed were prevented by the failure of 
the spinner retaining bolt which led to the loss of the 
stage. The cause of the failure was shown by analysis and 
tests will resume when replacement parts are delivered. 

A new data acquisition system, using PC-based control 
of VXI-bus instrumentation, was programmed and used 
successfully in the tests. HPVEE software was used with the 
VXI-bus hardware for Scanivalve control, acquisition of 
pressures, temperatures, strains and frequencies, and for 
the reduction of data. The system was found to be extremely 
flexible and well-suited to the researcher in the laboratory 
environment . 

The flow through the rotor alone was computed using the 
RVC3D code. The test results obtained agreed qualitatively 
with the predicted rotor performance and were also 
numerically consistent with them. This was a first step in 
the code validation process. The next step will be to 
measure the flow profile at the exit of the rotor and 
compare what is measured to the predicted flow field. The 
final step will be to measure the flow between the rotor 
blades and pressure distribution on the blade surfaces and 
coirpare these to predicted results. 
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with the experience obtained from the present study the 

following specific recommendations are made: 

• With the inlet tenperature fluctuations caused by the 
coastal location of the laboratory, increased use of 
averaging is required. First, a scan of temperatures and 
wheel speed should be conducted, then a pressure scan, 
then a second scan of tenperatures and wheel speed should 
be taken and averaged with the first one to define one 
data point. Three data points should then be taken and 
averaged to produce an averaged data point to plot on the 
performance maps . 

• The torque- flexure strain gauge bridge output should be 
readjusted to zero between each constant speed run. (The 
output should be zero with no torque applied. However, 
the zero tended to drift with increased operating 
temperatures and therefore needed to be adjusted) . 

• Install a strain gauge multiplexer for strain 
measurements. Software procedures would then make the 
previous recommendation unnecessary and would be the 
preferred solution. 

• In rebuilding the stage, use only a solid steel spinner 
retaining bolt, use a transparent case wall, and re-use 
the video-monitoring and recording system demonstrated in 
the present study. 
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APPENDIX A ANALYSIS OF THE SPINNER BOLT FAILURE 



A. OVERVIEW 

While operating at 80% speed, with the inlet throttled 
to produce near stall conditions, the bolt attaching the 
spinner to the rotor failed. The spinner detached and 
contacted the rotor, causing complete loss of the stage and 
casing. Photos of the failed bolt, the destroyed stage, and 
the damaged Lucite casing are shown in Figs. Al, A2, and A3 
respectively. While rxinning the con^ressor, a video camera 
was being used to monitor and film the operation. Using the 
video film (at 30 frames/sec) and pressure readings measured 
immediately before the event, an analysis was conducted of 
the failure. 

B . CONSTRUCTION 

The spinner used with the stage had been modified from 
the design used on the Vavra stage. The retaining bolt was 
required to be two inches longer; namely, 3.09 inches versus 
1.09 inches for the Vavra spinner retaining bolt. The bolt 
was a Socket Head Cap Screw #8-32, which was not 
manufactured (standard) with the new required length. In 
order to increase the length of the bolt, a 1-1/8 inches 
bolt was cut in half. A metal sleeve with 5/32 inches o.d. 
and 1/8 inches i.d. was placed over the two halves that had 
been turned down to 1/8 inches o.d. The metal sleeve was 
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then silver brazed at 1050'’ F to form the extension. After 

spin balancing the rotor and drive shaft, the spinner was 

attached using the bolt and the assembly rebalanced. This 
required removing material from the inner edge of the 

spinner. When balance was finally achieved, the bolt 
attaching the spinner was 'lightly torqued' . Only after the 
event was it found that the metal was cut from copper 
tubing . 

C. ANALYSIS 

In order to estimate the aerodynamic load on the 

spinner the pressure measurements recorded just before the 
failure occurred were used. The throttle had been closed a 
small fraction more when the failure occurred; therefore, 
the previously scanned pressure values are considered to be 
very close to those that were present . For a worst-case 
analysis, the pressure inside the spinner was assumed to be 
the stagnation pressure at the nose of the spinner. {This 
assumption is reasonable since the hole that the bolt passed 
through was not sealed and leakage in would have occurred. 
The shoulder at the rim of the spinner was also not sealed 
and the pressure behind the spinner would leak out and be 
less than this 'worst case' assumption). The pressure on the 
outer wall was taken to be the average of the static 
pressures taken at station one, another static pressure 
reading taken mid-length of the spinner, and the static 
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pressure at the nose of the spinner. All the static pressure 
measurements were taken at the casing wall. With this 
pressure distribution an axial force of 11.0 Ibf was 
calculated. This axial force would have caused a tensile 
stress of 1.6 ksi in the copper sleeve of the bolt. The 
stress in the copper sleeve due to the torque used to attach 
the bolt to the rotor was estimated to have been 10.3 ksi. 
The maximum steady total tensile stress in the copper sleeve 
was therefore 11.9 ksi. 

From the videotape of the event, it was observed that 
the spinner went upstream and disappeared from view within 2 
video frames, came back and 'flew' stably for several 
seconds, moving around the annulus, and then fell into the 
blading. Since the frame speed was known to have been 3 0 
frames/sec, an estimate of the acceleration of the spinner 
as it left the rotor was made by comparing the relative 
positions of the spinner over two frames. The estimated 
average acceleration was 225 ft/sec". The drag force on the 
spinner (downstream) was calculated to have been 2.52 Ibf. 
The net force required to accelerate the spinner upstream 
was calculated to have been 2.33 Ibf. The initial force on 
the spinner necessary to explain the video images was 
therefore 4.85 Ibf. This compares to the estimate made using 
pressure measurements that calculated a force of 11.0 Ibf. 
The image on the first frame following detachment was 
blurred as the spinner moved rapidly upstream, making an 
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accurate estimation of the instantaneous acceleration at the 
time of detachment impossible. Therefore, although the 
calculations do not agree closely, they do demonstrate that 
the forces calculated to be present were of the order of 
magnitude necessary to explain the observed behavior of the 
spinner if the copper sleeve failed in tension. 

The tensile strength of annealed copper is 32 ksi and 
the yield strength is 5 ksi. Copper has the property of 
strain hardening. Therefore, each time the sleeve was 
stressed into plastic deformation the tensile strength would 
have decreased. With a stress twice the yield strength 
present, failure of the bolt occurred when the tensile force 
exceeded the reduced tensile strength. 
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Fig. A1 View of the Failed Spinner 
Retaining Bolt 
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Fig. A2 View of Stage After Failure 
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Fig. A3 View of Lucite Casing After Failure 
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APPENDIX B RESULTS 



The reduced results of all test runs are presented 
below. Data points with no values present for horsepower 
were caused by problems with the data acquisition system. 
The point at the same state was re-measured until a complete 
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APPENDIX C SOFTWARE USER'S GUIDE 



A. HPVEE 

The procedure for operating the HPVEE data acquisition 
routine was as follows: 

1. The input array for the Scanivalve Control user's 

object was initialized. The element was the 

Scanivalve Position (1 through 4) and the next 48 
elements of the array indicated whether the 
corresponding port should be measured (0 - no skip, 
1 - yes measure) 

2. When the rig was stable and at the desired condition 
for measurement the START button was depressed. 
Desired quantities were input as prompted. 

3 . When the data acquisition process was finished the 
data were automatically written to the hard disk. 

4. The operating condition was changed and the process 
repeated from step 1. 

All programs used were saved on a 3.5 inch floppy and 
are in the possession of Dr. Shreeve. A graphical listing of 
the program, with redundant portions for scanivalve 
positions 2 through 4 removed, follows: 
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C:\BART\lcr_reduction_printvee 
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^ Formula ; sum(x) 

20 at: 07 Formula "Formula” 
^ Formula : A/B*C 
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Main Compressor calculations (M.10B): Detail View, Total Objects 45 
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I Si |jg£[| Compressor caicu iailonB 
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11 at: B6 Formula "Formula 
Formula: A*9/5 + 32 

12 at: A2 Start "Start" 



Error Message: "Remember, inches of mercury or 
at: D2 UserObject "Scanivalve Control Object" 
at: C5 AlphaNumeric "SV in volts" 
at; D4 Formula "Pressure ms Hg absolute" 
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PanelView "Pressure Measurement" (w=lOOO, h=802) 
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1997 12:23:32 




Error Improper Scanivalvo Chosen- 
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Sat 20/Sep/1997 12 23:32 

Main. P,esau,e Measurement. Scanivalve Control Object (M 197 4) Detail View. Total Objects 10 
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Main; Pressure Measurement: Scanivalve Control Object: SV Reset; Reset SV1 (M.197 4 29 1): Detail View, Total Objects; 7 
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C:\BARTMcr_reductionj)rint vee 

Sal 20/Sep/1997 12 23:32 

Main; Pressure MeasuremenI: Scanivalve Control Object SV 1 Data Acquisition: Step SV1 (M 1 97.4 45 1 3): Detail View. Total Objects 7 
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C:\BART\Icr_reducl«on_pnnl vee 
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Main: Pressure Display: Panel View 
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B. RVC3D 



All software used and results obtained using RVC3D and 
associated programs are stored in a directory 
/ jedi/d3/blgrossm on the NFS Cray. Dr. G. V. Hobson has 
access to the files. The procedure for operating RVC3D was 
as follows: 

1. Input the geometry provided as geometry.dat into 
TCGRID to obtain the grid file fort.l. geometry. out was 
produced and examined to ensure the grid was properly 
produced. The command used was "tcgrid.exec <geometry . dat 
>geometry . out" . 

2. Modified rvc.in to reflect the desired initial 
conditions to be analyzed [Ref. 3]. 

3. Srobmitted the job to the queue using rvc3d.nqs. 

4. After completion, a solution file q.sol and a 
result file rvc . out were output. 

5. If further iterations were desired, modified irvc.in 
such that iresti=l and renamed q.sol to q.res. 

6. Resubmitted the job to the queue. 

Ample assistance was available from cortputer center 
personnel. Examples of geometry.dat, rvc.in and rvc3d.nqs 
are included in the following pages . 
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Listing of rvc3d.nqs; 



#QSUB -q reg -IT 15280 -It 15280 -IM 9Mw -me 
set RUN = grid 

assign -F til -N ieee -a ${RUN}.ieee u:l 
set RUN = q 

assign -F til -N ieee -a ${RUN}.sol u:3 
assign -F til -N ieee -a ${RUN}.res u:2 



cd /d3/blgrossm/termproj 
rvc3d <rvc.in >rvc.out 



Listing of rvc.in: 

•sanger rotor* 

Sciill im=185 jm=40 km=49 itl=20 iil=85 Scend 

&:nl2 nstg=4 cfl=5 aviscl=0 avisc2 = 1.0 avisc4=.5 irs=l epi=0.60 epj=0.70 
epk=0.70 ivdt=l itmax=5100 &end 
&jil3 ibcin=l ibcex=3 ires = 10 icmt=50 iresti=l &end 
&jil4 om=-2.54 prat=1.1905 emxx=.46 alex=52 igeom=l fi^end 
&nl5 ilt=2 iltin=2 renr=6.651e6 tw=0 vispwr=.6667 cmutm=14 jedge=25 
kedge=18 srtip=0 xrle=-.2 xrte=0.3 dblh=0.001 dblt=0.02 &end 
SciilS &end 



Listing of geometry.dat: 

itnaml merid=3 im=185 jm=40 km=49 itl=20 icap=16 k2d=ll Scend 
&Jiaiti2 nle=15 nte=10 dsle=3.5e-4 dste=2.5e-4 dshub=1.5e-5 dstip=1.5e-5 
dswte=3.5e-4 dswex=.010 dsthr=.39 dsmin=1.2e-5 dsmax=.002 
dsin=.010 clear=. 00000 rcorn=.008 Scend 
&nam3 iterm=100 idbg=0 0000000 aabb= . 6 ccdd=.2 &end 
&nam4 zbc=-0.100 -0.040 0.2000 -0.1000 -0.0400 0.2000 

rbc= 0.220 0.220 0.2800 0.4560 0.4560 0.4560 &end 

‘NPGS TRANSONIC DESIGN -GENER. FROM V8FBND13 INCDP-BASE=DEC13 INMER* 

32 14 

-.833 -.667 -.5833 -.50 -.4583 -.425 -.40833 -.3917 
-.3333 -.25 -.16666 -.0833 0. .0833 .1 .1167 

.125 .133 .14167 .15 .1583 .16666 .175 .183 
.19166 .208 .25 .3333 .41666 .5 .667 .8333 

.00083 .00083 .00083 .00083 .001667 .00625 .01083 .01667 
.04467 .08937 .13405 .17873 .2234 .2681 .277 .286 
.2904 .29333 .295583 .297667 .29933 .30033 .30166 .3025 
.303333 .304166 .306083 .306083 .306083 .306083 .306083 .306083 

-.8333 -.7637 -.666666 -.5 -.33333 -.2083 -.1239 -.04167 
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.0529 .2167 .3083 .5 .66666 .83333 



.4583 .4583 .4583 .4583 .4583 .4583 .4583 .4583 
.4583 .4583 .4583 .4583 .4583 .4583 



11 23 22 
0.000000 
0.025568 
0.031878 
0.081719 
0.089506 
0.001360 
0.026281 
0.032350 
0.079987 
0.087529 
0.002013 
0.026682 
0.032644 
0.078992 
0.086358 
0.002617 
0.027000 
0.032848 
0.077748 
0.084894 
0.003875 
0.027672 
0.033327 
0.075885 
0.082643 
0.007383 
0.029584 
0.034782 
0.072754 
0.078711 
0.010925 
0.031430 
0.036196 
0.070365 
0.075575 
0.012591 
0.032161 
0.036707 
0.069210 
0.074074 
0.014092 
0.032945 
0.037327 
0.068699 
0.073344 
0.014068 
0.032800 
0.037163 
0.068472 
0.073071 



0.000514 

0.038395 

0.097415 

0.001751 

0.038563 

0.095222 

0.002353 

0.038731 

0.093887 

0.002913 

0.038800 

0.092212 

0.004123 

0.039064 

0.089584 

0.007577 

0.040025 

0.084862 

0.011094 

0.040990 

0.081001 

0.012749 

0.041277 

0.079161 

0.014243 

0.041734 

0.078207 

0.014213 

0.041553 

0.077861 



0.002301 

0.045135 

0.104086 

0.003532 

0.044956 

0.101740 

0.004135 

0.044945 

0.100278 

0.004692 

0.044857 

0.098434 

0.005878 

0.044881 

0.095503 

0.009239 

0.045312 

0.090132 

0.012640 

0.045811 

0.085687 

0.014226 

0.045871 

0.083566 

0.015665 

0.046166 

0.082428 

0.015624 

0.045971 

0.082024 



0.005027 

0.052088 

0.109470 

0.006238 

0.051567 

0.107019 

0.006837 

0.051351 

0.105463 

0.007384 

0.051036 

0.103490 

0.008529 

0.050778 

0.100323 

0.011743 

0.050643 

0.094442 

0.014965 

0.050659 

0.089541 

0.016447 

0.050490 

0.087197 

0.017803 

0.050623 

0.085913 

0.017747 

0 . 050*417 

0.085465 



0.008746 

0.059237 

0.113531 

0.009908 

0.058387 

0.111013 

0.010492 

0.057968 

0.109392 

0.011017 

0.057417 

0.107325 

0.012096 

0.056765 

0.103986 

0.015099 

0.056018 

0.097727 

0.018075 

0.055534 

0.092492 

0.019417 

0.055133 

0.089983 

0.020663 

0.055105 

0.088589 

0.020588 

0.054890 

0.088111 



0.013525 

0.066568 

0.116248 

0.014595 

0.065405 

0.113692 

0.015144 

0.064787 

0.112029 

0.015627 

0.064001 

0.109902 

0.016608 

0.062939 

0.106450 

0.019322 

0.061437 

0.099942 

0.021980 

0.060436 

0.094488 

0.023145 

0.059800 

0.091869 

0.024255 

0.059611 

0.090402 

0.024157 

0.059390 

0.089904 



0.019448 

0.074067 

0.116636 

0.020362 

0.072609 

0.114041 

0.020848 

0.071798 

0.112356 

0.021259 

0.070780 

0.110207 

0.022099 

0.069314 

0.106724 

0.024430 

0.066995 

0.100171 

0.026692 

0.065366 

0.094690 

0.027641 

0.064492 

0.092059 

0.028587 

0.064143 

0.090584 

0.028465 

0.063918 

0.090081 
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0.014718 
Q. 032912 
0.037239 
0.067997 
0.072518 

0.222890 

0.236541 

0.239910 

0.266519 

0.270677 

0.251549 

0.262257 

0.264864 

0.285332 

0.288573 

0.267527 

0.276842 

0.279093 

0.296593 

0.299375 

0.284586 

0.292532 

0.294438 

0.309070 

0.311399 

0.309498 

0.315614 

0.317067 

0.328005 

0.329743 

0.357485 

0.360641 

0.361380 

0.366779 

0.367626 

0.398682 

0.400090 

0.400417 

0.402763 

0.403120 

0.418757 

0.419585 

0.419778 

0.421154 

0.421360 

0.428363 

0.428980 

0.429124 

0.430151 

0.430303 

0.441651 

0.442033 

0.442122 

0.442760 

0.442854 



0.018291 

0.050230 

0.084538 

0.225574 

0.250699 



0.021057 

0.054624 

0.087114 

0.227560 

0.254516 



0.014854 

0.041537 

0.077153 

0.223164 

0.243389 

0.274899 

0.251717 

0.267534 

0.291878 

0.267655 

0.281391 

0.302218 

0.284683 

0.296378 

0.313784 

0.309561 

0.318542 

0.331527 

0.357513 

0.362126 

0.368500 

0.398694 

0.400746 

0.403493 

0.418763 

0.419971 

0.421576 

0.428368 

0.429268 

0.430462 

0.441654 

0.442211 

0.442951 



0.016226 

0.045868 

0.081192 

0.224119 

0.246988 

0.278461 

0.252483 

0.270281 

0.294678 

0.268328 

0.283738 

0.304631 

0.285263 

0.298351 

0.315811 

0.310012 

0.320037 

0.333048 

0.357749 

0.362877 

0.369249 

0.398800 

0.401077 

0.403815 

0.418826 

0.420166 

0.421762 

0.428414 

0.429413 

0.430600 

0.441683 

0.442301 

0.443036 



0.281335 

0.253645 



0.296946 

0.269348 



0.306589 

0.286140 



0.317459 

0.310694 



0.334287 

0.358105 



0.369862 

0.398959 



0.404079 

0.418920 



0.421916 

0.428484 



0.430714 

0.441726 



0.283503 

0.255222 



0.298663 

0.270728 



0.308072 

0.287324 



0.318709 

0.311611 



0.335228 

0.358582 



0.370329 

0.399173 



0.404282 

0.419046 



0.422034 

0.428578 



0.430802 

0.441784 



0.024535 

0.059050 

0.088863 

0.230111 

0.258430 

0.284954 

0.257236 

0.279067 

0.299814 

0.272485 

0.291230 

0.309068 

0.288826 

0.304590 

0.319549 

0.312770 

0.324678 

0.335861 

0.359182 

0.365170 

0.370644 

0.399441 

0.402081 

0.404419 

0.419204 

0.420756 

0.422114 

0.428696 

0.429853 

0.430861 

0.441857 

0.442575 

0.443197 



0.028737 

0.063508 

0.089032 

0.233273 

0.262434 

0.285162 

0.259714 

0.282162 

0.299964 

0.274639 

0.293877 

0.309192 

0.290661 

0.306799 

0.319648 

0.314181 

0.326317 

0.335932 

0.359909 

0.365960 

0.370677 

0.399764 

0.402419 

0.404433 

0.419394 

0.420955 

0.422122 

0.428837 

0.430002 

0.430867 

0.441945 

0.442667 

0.443201 



0.442392 0.442483 
0.443106 0.443160 



0.273121 0.276052 



0.286156 0.288655 



0.300365 0.302444 



0.321552 0.323091 



0.363635 0.364399 



0.401410 0.401744 



0.420362 0.420558 



0.429559 0.429706 
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0.455686 

0.455812 

0.455841 

0.456054 

0.456085 

0.000000 

0.130240 

0.155860 

0.288078 

0.300126 

0.003128 

0.121816 

0.146442 

0.287196 

0.302128 

0.006040 

0.119448 

0.143479 

0.287186 

0.303358 

0.007812 

0.116621 

0.140101 

0.286148 

0.303448 

0.009343 

0.113166 

0.136014 

0.284567 

0.303219 

0.012457 

0.109937 

0.131892 

0.281254 

0.301299 

0.015765 

0.109331 

0.130593 

0.277785 

0.298246 

0.017368 

0.109313 

0.130251 

0.275683 

0.296131 

0.017681 

0.109172 

0.130021 

0.274953 

0.295433 

0.020481 

0.110394 

0.130889 

0.273383 

0.293579 



0.455711 

0.455931 

0.456168 

0.029582 

0.222896 



0.455730 

0.455961 

0.456186 

0.050062 

0.241772 



0.455687 

0.455871 

0.456117 

0.003130 

0.179937 

0.310517 

0.005248 

0.170204 

0.315751 

0.007780 

0.166850 

0.318393 

0.009253 

0.163087 

0.319781 

0.010493 

0.158540 

0.321104 

0.013335 

0.153720 

0.320851 

0.016553 

0.151800 

0.318337 

0.018125 

0.151148 

0.316262 

0.018424 

0.150833 

0.315629 

0.021191 

0.151350 

0.313542 



0.455696 

0.455901 

0.456145 

0.013824 

0.202284 

0.317930 

0.014776 

0.192927 

0.326106 

0.016799 

0.189543 

0.330054 

0.017833 

0.185579 

0.332650 

0.018600 

0.180744 

0.335414 

0.020849 

0.175419 

0.336758 

0.023727 

0.172950 

0.334786 

0.025166 

0.172006 

0.332787 

0.025428 

0.171610 

0.332229 

0.028075 

0.171776 

0.329959 



0.323057 

0.028894 



0.333741 

0.030195 



0.338813 

0.030604 



0.342455 

0.030696 



0.346469 

0.032095 



0.349226 
0. 034476 



0.347746 

0.035718 



0.345836 

0.035926 



0.345354 

0.038395 



0.326440 

0.047394 



0.339090 

0.047807 



0.345050 

0.047445 



0.349523 

0.046703 



0.354528 

0.047041 



0.358424 

0.048786 



0.357349 

0.049774 



0.355522 

0.049911 



0.355104 

0.052140 



0.455754 

0.455992 

0.456198 

0.074802 

0,258919 

0.328478 

0.069991 

0.253398 

0.342478 

0.069419 

0.251457 

0.349050 

0.068195 

0.248690 

0.354099 

0.066514 

0.245014 

0.359790 

0.065645 

0.239739 

0.364482 

0.066639 

0.236060 

0.363693 

0.067318 

0.234336 

0.361930 

0.067370 

0.233725 

0.361559 

0.069302 

0.232846 

0.358983 



0.455783 

0.456023 

0.456199 

0.103219 

0.274350 

0.328755 

0.096331 

0.270953 

0.342908 

0.094760 

0.269884 

0.349537 

0.092652 

0.267892 

0.354630 

0.089999 

0.265160 

0.360368 

0.087855 

0.260729 

0.365102 

0.088013 

0.256983 

0.364329 

0.088335 

0.255031 

0.362568 

0.088289 

0.254357 

0.362199 

0.089865 

0.253132 

0.359607 



0.192168 0.212524 
0.342944 0.352593 



0.214376 0.234537 



0.211273 0.231916 



0.207502 0.228555 



0.202624 0.224143 



0.196989 0.218430 



0.194043 0.215080 



0.192823 0.213600 



0.192351 0.213056 
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0.022265 

0.111137 

0.131402 

0.272292 

0.292279 

0.012140 

0.046930 

0.052625 

0.039421 

0.032344 

0.008415 

0.038449 

0.043927 

0.037468 

0.030935 

0.006996 

0.034788 

0.040080 

0.036547 

0.030278 

0.005881 

0.031279 

0.036310 

0.036021 

0.030002 

0.004812 

0.027160 

0.031787 

0.035758 
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